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METHOD FOR GEOLOGIC MODELING THROUGH HYDRODYNAMICS- 
BASED GRTODING (HYDRO-GRIDS) 

This application claims the benefit of U.S. Provisional Application No. 
5 60/584,617 filed on July 1, 2004. 

FIELD OF THE INVENTION 

[0001] This invention relates generally to the fields of geophysical prospecting 
10 and reservoir characterization. More specifically, the invention is a method for 
creating geologic models permitting enhanced representation of the properties of a 
subsurface volume, such as a reservoir or basin. 

* 

BACKGROUND OF THE INVENTION 

[0002] A geologic model is a computer-based representation of a subsurface 

15 earth volume, such as a petroleum reservoir or a depositional basin. In the oil and 
gas industry, geologic models are commonly used in activities such as determining 
the locations of wells, estimating hydrocarbon reserves, or planning reservoir- 
development strategies. Geologic models are also input to fluid flow simulations to 
test development and production scenarios in order to optimize return on investment 
20 A key parameter in flow simulation is the spatial distribution of permeability, which 
together with the properties of the hydrocarbons and other fluid found in the 
subsurface reservoir, determines the producibility of the reservoir. Typically, the 
input data for geologic models includes data obtained through seismic prospecting 
techniques. 

25 [0003] Seismic prospecting techniques are commonly used to aid in the search 

for and evaluation of subterranean hydrocarbon reservoirs. A seismic prospecting 
operation consists of three separate stages: data acquisition, data processing, and data 
interpretation. The success of the operation depends on satisfactory completion of all 
three stages. 

30 [0004] In the data acquisition stage, a seismic source is used to generate an 

acoustic signal that propagates into the earth and is at least partially reflected by 
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subsurface seismic reflectors. The reflected signals are detected and recorded by an 
array of seismic receivers located at or near the surface of the earth, in an overlying 
body of water, or at known depths in boreholes. 

[0005] During the data processing stage, the recorded seismic signals are 

5 refined and enhanced using a variety of procedures that depend on the nature of the 
geologic structure being investigated and on the characteristics of the raw data. In 
general, the purpose of the data processing stage is to produce an image of the 
subsurface from the recorded seismic data for use during the data interpretation stage. 

[0006] The purpose of the data interpretation stage is to determine information 

10 about the subsurface geology of the earth from the processed seismic data. The 
results of the data interpretation stage may be used to determine the general geologic 
structure of a subsurface region, or to locate potential hydrocarbon reservoirs, or to 
guide the development of an already discovered reservoir. 

[0007] At present, the conclusions, which can be made after the data 

15 interpretation stage, are generally limited to broad descriptions of the size and general 
nature of subsurface structures due to the limited resolution of seismic data. The 
descriptions may for example provide an indication of the total volume of 
hydrocarbons, which might be retained in such structures. However, present 
• technology does not allow the analyst to determine production rates from the 
20 subsurface formations if hydrocarbons are discovered. In addition, when an 
exploration well has been drilled, present technology does not allow an analyst to be 
able to accurately characterize the nature of the subsurface internal geometry in 
locations other than in the most immediate region of any such well. In particular, 
reservoir permeability and continuity are not well characterized by present 
25 technology. 

[0008] The hydrocarbon volume and rate of production depend on a variety 

of factors, including fluid properties, reservoir net-to-gross, porosity, permeability, 
spatial variability in grain size distribution, and connectivity. Reservoir connectivity, 
a measure of the communication (or lack thereof) between points within the reservoir, 
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is a strong function of the reservoir internal geometry and is commonly a primary 
factor controlling hydrocarbon production efficiency and ultimate recovery. There is 
a need to predict the detailed internal geometry of subsurface reservoirs using 
geologic data, such as seismic data, and without having to drill many exploration and 
5 delineation wells. Such a capability would facilitate estimation of hydrocarbon 
volume in place and production rates early in the hydrocarbon exploration and 
development process. 

[0009] Another step in geologic modeling is "gridding Gridding is the 

division of the subsurface region into cells, within which the rock properties are 

10 regarded as uniform. The ability to accurately model reservoir internal geometry and 
connectivity is largely dependent on the vertical grids used in the modeling processes 
due to the nature of stratification and vertical heterogeneity in the reservoir. As a 
result, the gridding methods in geologic modeling focus mainly on the vertical 
variation of attributes of a geologic model. In map view, uniform rectangular grids 

15 are commonly used for relatively small variations of attributes in the lateral 
directions. 

[0010] There are four gridding methods commonly used in commercially 

available geologic modeling tools. Examples of commercially available geologic 
modeling tools include: SGM (Landmark Graphics Corporation, Stratamodel 

20 Geocellular Modeling (SGM), 1989-2003), RMS (Roxar ASA, Irap RMS (Reservoir 
Modeling Systems) , 1993-2003), gOcad (Earth Decision Sciences Corporation, 
gOcad, 1989-2003), and Petrel (Technoguide (A Schlumberger Product Group), 
Petrel*" Worfylow Tools, 1996-2003). The four gridding methods used in the 
petroleum industry are commonly named proportional, onlap, truncation, and 

25 reference. The name of each method can vary from one tool to another. For example, 
"reference-grid" in RMS is the same as "depositional-grid" in SGM. The choice of the 
most appropriate gridding method to use when building a geologic model depends on 
the geologic modeler's experience and personal judgement of the model's particular 
geologic setting. Gridding style, however, was created to mimic natural patterns of 

30 erosion and deposition observed in nature at the seismic scale. 
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[0011] The gridding methods used in the most commercially available 

geologic modeling tools (such as, SGM, RMS, gOcad, and Petrel) are based on 
methods disclosed in two of Swanson's patents. The patents are U.S. Patent Nos. 
4,821,164 and 4,991,095. 

5 [0012] Figures 1 (a) is a cross-section illustration of stratification of sand 

bodies 1, 2, 3, 4 and 5 and low permeability layers 6. Figures 1(b), 1(c), 1(d), and 
1(e) illustrate sand bodies 1, 2, 3, 4, and 5 and low permeability layers 6 from figure 
1(a) using conventional grids. 

[0013] As shown in figure 1(b), "Proportional-grid" assumes that the reservoir 

10 sand bodies 1, 2, 3, 4, and 5 within the given sequence are parallel to the top 7 and 
basal surfaces 9 of the sequence 11. Therefore, the sequence is subdivided into the 
cells 13 proportional to the thickness of the sequence 11. 

[0014] As shown in figure 1(c), "Onlap-grid" assumes that the sand bodies 1, 

2, 3, 4, and 5 are parallel to the top surface 7 of the sequence 11. The sequence 11 is 
15 gridded with the constant thickness cells 13 that are parallel to the top surface 7 and 
may be truncated 15 by the basal surface 9 of the sequence 11. 

[0015] As shown in figure 1(d), "Truncation-grid" assumes that the sand 

bodies 1, 2, 3, 4, and 5 are parallel to the basal surface 9 of the sequence 11. The 
sequence 11 is gridded with the constant thickness cells 13 that are parallel to the 
20 basal surface 9 and may be truncated 15 by the top surface 7 of the sequence 11. 

[0016] As shown in figure 1(e), "Reference-grid" assumes that the sand 

bodies 1, 2, 3, 4, and 5 are parallel to a given surface (or reference surface) 17, such 
that the constant thickness cells 13 are constructed parallel to the reference surface 17 
and may be truncated 15 by the top 7 and basal surfaces 9 of the sequence 11. The 
25 reference surface 17 is usually some sort of geologic datum. In this figure, the 
reference surface is the bottom of sand body 2. 

[0017] The problem is that none of these gridding methods represent the true 

internal geometry accurately. To illustrate this problem, there are five disconnected 
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sand bodies 1, 2, 3, 4, and 5 within the sequence 11 of figure 1(a) and the above 
gridding methods create significantly distorted representations as shown in figures 
1(b), 1(c), 1(d), and 1(e) when compared to figure 1(a). Sand body 1 is incorrectly 
connected with sand body 2 in the proportional-grids and onlap-grids as illustrated in 
5 figures 1(b) and 1(c) respectively. As shown in figure 1(c), sand body 2 is incorrectly 
connected with sand body 3 in the onlap-grids. Furthermore, sand body 4 is 
incorrectly connected with sand body 5 in the truncation-grids of figure 1(d). 

[0018] Some sand bodies are incorrectly disconnected within themselves. For 

example, sand body 4 is disconnected in the proportional-grids of figure 1(b) and sand 

10 body 2 is disconnected in the onlap-grids of figure 1(c). Furthermore, some sand 
bodies are erroneously disconnected from the fluid flow since some cells in these sand 
bodies only contact their neighbor cells at the corners, without face connections. For 
example, sand body 2 is disconnected in the proportional-grids of figure 1(b) and sand 
body 4 is disconnected in all of the four gridding methods as shown in figures 1(b), 

15 1(c), 1(d), and 1(e). Misalignment between the sand bodies and grids is the direct 
cause of this distortion problem. 

[0019] The grids constructed by the above-mentioned prior art methods do not 

properly align themselves with the reservoir internal geometry. The alignment 
between the grids and the reservoir internal geometry is important for two reasons. 

20 First, any lack of alignment distorts the reservoir connectivity in the geologic model 
subsequently built Second, when the geologic model is scaled up for reservoir 
performance simulation with coarser grids, the coarser grids will magnify the 
misrepresentations and distortions that exist in the geologic model. Consequently, the 
accuracy of any reservoir performance prediction will decrease, which will hamper 

25 reservoir management. 

[0020] A reliable reservoir performance simulation depends strongly on the 

ability of a geologic model to accurately characterize the spatial distribution of 
permeability. Most important are the permeability extremes, such as high- 
permeability sand bodies and no-flow shale barriers, because the permeability 
30 extremes typically control the oil, gas, and water flows. The ability to model 
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permeability is currently severely limited by the commercially available geologic 
model gridding methods. One reason that the existing geologic modeling tools have 
to adopt arbitrary lateral grid surfaces is that the actual strata geometry (or layering) 
of the sand objects is not known a priori. As a consequence, the reservoir 
5 connectivity is typically misrepresented in geologic models and the resulting 
performance simulation models. Accordingly, there is a need for a method that 
generates grids that follow the geometry of sand bodies and shale barriers. 
Preferably, these grids would allow the reservoir connectivity to be accurately 
characterized in a geologic model and preserved during upscaling for reservoir 
10 performance simulation. The current invention addresses the need by generating the 
lateral grid surfaces naturally, following the complex depositional and erosional 
processes observed in nature. 

SUMMARY OF THE INVENTION 

[0021] Under the present techniques, various embodiments of a hydrodynamics- 
15 based gridding method, called 'Hydro-Grids", may generate grid surfaces conforming 
to time surfaces during the process of sedimentation, which create a subsurface 
sedimentary volume. The location of these time surfaces is typically determined by 
sedimentary process simulation and interpretation of geologic data. Since these 
surfaces can precisely follow the geometry of sand bodies and shale barriers, geologic 
20 models conforming to these surfaces may accurately characterize sedimentary 
connectivity and may preserve their accuracy during upscaling for reservoir 
performance simulation. 

[0022] In one embodiment, a method comprises obtaining geologic data, creating 
horizontal grid lines wherein at least one horizontal grid line represents deposit time 

25 lines based on the deposition rates from the geologic data, creating vertical grid lines, 
and modeling the sedimentary properties of the sedimentary basin or reservoir from 
the grids of geologic data created by the vertical and horizontal grid lines. The 
geologic data used includes but is not limited to seismic data, outcrop studies, well log 
data, core data, forward and inverse modeling data, correlation of seismic data to 

30 sedimentary data. The modeled sedimentary properties of the reservoir include, but 
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are not limited to grain size distribution, connectivity, net-to-gross, porosity, 
permeability, pore pressure and any combination thereof. 

[0023] A second embodiment for geologic modeling of a sedimentary basin or 

reservoir is also disclosed. This embodiment comprises obtaining geologic data, 
5 creating at least one horizontal grid line based on an estimated time from the 
deposition and erosion in the geologic data, determining stratigraphic hierarchy levels 
of the horizontal grid lines, merging neighboring horizontal grid lines, upscaling the 
merged layers based on the stratigraphic hierarchy, estimating sedimentary properties 
from the upscaled merged layers. The hierarchy levels assigned to the horizontal grid 
10 lines are selected from methods including but not limited to changes in the inlet 
conditions, changes in the simulation conditions, variations in the vertical distribution 
of the composition of the simulated deposit and any combination thereof. 

■ 

[0024] A method for geologic modeling of a subsurface volume, such as a 
geologic basin or reservoir, is disclosed. In one embodiment the method comprises 

15 obtaining geologic data, creating vertical grid surfaces to divide the subsurface 
volume into "columns", creating lateral grid surfaces to correspond with surfaces of 
constant geologic time during the creation of the geologic structures in the subsurface 
volume, and representing the geologic properties of the subsurface volume within the 
grids cells bounded by the vertical and lateral grid surfaces. The geologic data used 

20 may include, but is not limited to, seismic data, outcrop studies, well log data, core 
data, data from numerical modeling of geological processes, stratigraphic surfaces 
interpreted from seismic data. The modeled geologic properties of the reservoir may 
include, but are not limited to, grain size distribution, connectivity, net-to-gross, 
porosity, permeability, pore pressure and any combination thereof. 

25 [0025] A second embodiment for geologic modeling of a subsurface volume is 

also disclosed. This embodiment comprises obtaining geologic data, creating vertical 
grid surfaces to divide the subsurface volume into "columns", creating lateral grid 
surfaces to correspond with surfaces of constant geologic time during the creation of 
the geologic structures in the subsurface volume, representing the geologic properties 

30 of the subsurface volume within the grid cells bounded by the vertical and lateral grid 
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surfaces, merging selected grid cells to reduce the number of cells in the model, and 
establishing new values for the geologic properties in the merged grid cells. 

BRIEF DESCRIPTION OF THE DRAWINGS 

5 [0026] Figure 1(a) is a cross-section illustration of original sand bodies; 

[0027] Figure 1(b) is a cross-section illustration of proportional-grids 

representation of the original sand bodies in figure 1(a); 

[0028] Figure 1(c) is a cross-section illustration of onlap-grids representation 

of the original sand bodies in figure 1(a); 

10 [0029] Figure 1(d) is a cross-section illustration of truncation-grids 

representation of the original sand bodies in figure 1(a); 

■ 

[0030] Figure 1(e) is a cross-section illustration of reference-grids 

representation of the original sand bodies in figure 1(a); 

[0031] Figure 2 is a flow chart of a first embodiment of the invention; 

15 [0032] Figure 3 is a top view of 4671 unstructured-cell deposit that would 

need 24806 structured grid cells for the same resolution; 

[0033] Figure 4(a) is a cross-section illustration of original sand bodies similar 

to figure 1(a); 

[0034] Figure 4(b) is a cross-section illustration of the sand bodies from figure 

20 3(a) in Hydro-Grids; 

[0035] Figure 5 is a flow chart of a second embodiment of the invention; 

[0036] Figure 6(a) is a hierarchy arrangement structure in a Hydro-Grids 

embodiment; 

[0037] Figure 6(b) is a hierarchy arrangement of geologic layers; 
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[0038] Figure 7 is a Hydro-Grids cross-section of a complex avulsive sand 

body composed of several sand body lobes and channels; 

[0039] Figure 8(a) is a cross-section illustration of sand body locations and 

sizes, and example deposit time surfaces populated by a sedimentary process network 
5 simulator; 

[0040] Figure 8(b) is a cross-section illustration of deposit time surfaces on 

fan delta populated by a sedimentary process network simulator. 

■ 

DETAILED DESCRIPTION 

[0041] In the following detailed description, the invention will be described in 

10 connection with its preferred embodiments. However, to the extent that the following 
description is specific to a particular embodiment or a particular use of the invention, 
this is intended to be illustrative only. Accordingly, the invention is not limited to the 
specific embodiments described below, but rather, the invention includes all 
alternatives, modifications, and equivalents falling within the true scope of the 
15 appended claims. 

[0042] Hydro-Grids are grids created by sedimentary process simulation. The 

Hydro-Grids lateral or horizontal surfaces, which may include and be represented by 
grid lines, are time surfaces recorded during the simulation of sand 
deposition/erosion. The grid surfaces follow the sand body shapes and their internal 
20 architecture (structure) naturally without any assumptions used by the existing 
geologic model gridding methods. Therefore, the sand body shapes and their internal 
structure can be precisely maintained for reservoir performance simulation. 

[0043] As shown in figure 2, one embodiment has four steps. In reference to 

numeral 21 of figure 2, geologic data is obtained. Now referring to numeral 22, 
25 vertical grid surfaces (lines) are created to divide the subsurface volume into 
"columns" of limited lateral extent. Now referring to numeral 23, horizontal or lateral 
grid surfaces (lines), which may be substantially or relatively horizontal, are created 
to correspond with surfaces of geologic time associated with the original deposition of 
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the sediments that fill the subsurface volume. Together, the vertical and lateral grid 
surfaces (lines) divide the subsurface volume into grid cells. Now referring to 
numeral 24, the geologic properties of the subsurface volume are represented by 
assigning to each grid cell a particular value for each geologic property. The geologic 
5 properties include but are not limited to grain size distribution, connectivity, net-to-, 
gross, porosity, permeability and pore pressure. 

[0044] Now referring to numeral 21 in figure 2, the geologic data is typically 

seismic data, but can also include core data, well logs, outcrop studies, time dating of 
deposits, numerical modeling results, geologic models, and interpreted stratigraphic 

10 surfaces. This data will typically include information that can aid in choosing the 
location of the surfaces, which may divide the subsurface volume into cells. 
Typically, this includes data on the lateral rates of change of geologic properties, 
which indicate how far apart vertical surfaces can be placed without losing accuracy. 
This also includes data that can aid in establishing boundary conditions for the 

15 numerical simulation of sedimentation which is part of the preferred embodiment for 
detenmning the spacing between the lateral surfaces in the step shown by numeral 23 
of figure 2. 

[0045] Now referring to numeral 22 in figure 2, vertical grid surfaces (lines) 

are selected to divide the subsurface volume into "columns". This can be done for 

20 both structured and unstructured gridding. In top or map view, for structured 
gridding, the vertical grid surfaces (lines) would divide the subsurface volume into 
uniform rectangular or square cells, while in unstructured gridding, the vertical grid 
surfaces (lines) could divide the volume into cells of varying shapes and sizes. The 
use of unstructured grids provides two advantages over structured grids. First, fine- 

25 scale cells can be constructed where the change in geologic properties is rapid, while 
coarse-scale cells can be used where the properties change slowly. As a result, the 
number of cells required to represent the subsurface volume with a particular level of 
accuracy should be significantly smaller with unstructured grid cells than with 
structured grid cells. For example, figure 3 shows a map view of 4,671 unstructured- 

30 cell deposit that would need 24,806 structured grid cells 33 for the same accuracy. 
Second, due to increased flexibility, unstructured grids capture deposit shapes very 
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accurately, while the structured grids may distort the deposit shapes with the same 
number of cells. These shapes include the delta fan deposit 35 and the channel 37 
cutting through the delta fan deposit 35. 

[0046] Now referring to numeral 23 in figure 2, the determination of lateral 

5 grid surfaces (lines) is typically accomplished with the aid of a numerical simulation 
of sedimentary processes. For example, SEDSIM, a product of CSIRO Petroleum 
Resources, is a suitable sedimentary process simulator that is commercially available 
and known to persons skilled in the art. The preferred embodiment uses the 
numerical simulation technique described in U.S. Provisional Patent Application No. 
10 60/608,707. 

[0047] Determining the lateral grid surfaces (lines) typically includes various 

steps. First, a basal surface for a subvolume within the subsurface volume is 
determined from seismic data. This is accomplished by a skilled seismic interpreter 
according to methods familiar to persons of ordinary skill in the ait Second, the basal 

15 surface is adjusted to correct its dip and orientation to the original position of the 
surface at the time of deposition. Depending on the nature and degree of post- 
depositional change, this may include a rotation of the surface to restore the 
paleoslope and unbending the surface to remove the effect of subsidence. The size 
and degree of these adjustments is determined according to the judgement of skilled 

20 interpreters as will be understood by persons of ordinary skill in the art. Computer 
simulation of subsidence and other post-depositional processes may assist in this 
determination (e.g., Bezemer, T.D., Kooi, H., and Cloetingh: "Numerical Modeling of 
Fault-Related Sedimentation," SEPM Special Publications No. 62, Numerical 
Experiments in Stratigraphy: Recent Advances in Stratigraphic and Sedimentologic 

25 Computer Simulations, SEPM (Society for Sedimentary Geology) 1999). Third, a 
numerical simulation of the sedimentary processes, which may fill the subsurface 
volume, is performed over the adjusted basal surface. Such simulations run forward 
in geologic time, computing the amount and nature of sediment deposited at each 
location at each time step. The top surfaces of the simulated deposit are recorded at 

30 regular time steps. The number of layers in the geologic model may be controlled by 
changing the time step size. Alternatively, the top surfaces may be recorded only at 
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time steps when a significant change in deposition rate or deposited sediment 
properties occurs. These steps may be iterated to adjust the boundary conditions of 
the simulation until good correspondence is achieved between the simulated deposit 
and the geologic data. Finally, if transformations were applied to the basal surface 
5 prior to simulation, the reverse of these transformations are applied to the basal 
surface and all computed top surfaces, so that the surfaces are moved to their likely 
present-day position. 

[0048] If, during the simulation, the flow is near steady state and 

predominantly depositional, the deposit composition is vertically homogeneous, and 
10 lateral grid surfaces may be separated by an extended period of simulation time. 
However, if the flow is highly transient with significant avulsions, the lateral grid 
surfaces may be defined at much smaller time intervals. This non-uniform-layering 
method ensures that the rock properties are as homogeneous as possible within a 
layer, which effectively characterizes reservoir connectivity. 

15 [0049] During and after the simulation for a particular subvolume, there are 

several adjustments that are typically made to the lateral grid surfaces (lines). First, 
the vertical coordinates of a lateral grid surface (line) may be adjusted to account for 
any removal of sediment by erosion after the original sediment in the layer is 
deposited. Otherwise the lateral grid surfaces (lines) may cross each other, which is 

20 not permitted. Second, it may be found that certain lateral grid surfaces (lines) 
created in simulation may not match in location the corresponding actual surfaces 
observed in seismic data, well logs, or other forms of geologic data as described in 
numeral 21 of figure 2. In this case, the lateral grid surfaces (lines) may be altered 
manually, or the boundary conditions for the simulation may be adjusted until a more 

25 favorable correspondence between the simulated and observed surfaces is obtained. 
The adjustment of boundary conditions is the preferred method when possible. 
However, in practice, this can be time consuming and may not produce an adequate 
correspondence, in which case manual adjustment may be made. 

[0050] Now referring to numeral 24 in figure 2, the geologic properties of the 

30 subsurface volume are represented by assigning values for these properties to each of 
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the grid cells created by the intersections of the vertical and lateral grid surfaces 
(lines). Typically, these values are determined by the same numerical simulation, 
which established the location of the lateral grid surfaces (lines). 

[0051] An example of Hydro-Grids is shown in figure 4. Figure 4(a) shows a 

5 cross-section of a simulated subsurface volume. In figure 4(b), the Hydro-Grids 
gridding method preserves the shapes of the five sand bodies 1, 2, 3, 4, and 5 between 
the low permeability layers 6. Because the Hydro-Grids surfaces (lines) follow the 
sand body shapes in figure 4(b), the connectivity between and within the original sand 
bodies in figure 4(a) is accurately represented. 

10 [0052] As shown in figure 5, one embodiment is now disclosed. Now 

referring to numeral 41 of figure 5, geological data is obtained. Now referring to 
numeral 42 of figure 5, vertical grid surfaces (lines) are created to divide the 
subsurface volume into "columns" of limited lateral extent. Now referring to numeral 
43 of figure 5, lateral grid surfaces (lines), which may be substantially or relatively 

■ 

15 horizontal, are created which correspond with surfaces of geologic time associated 
with the original deposition of the sediments that fill the subsurface volume. Together, 
the vertical and lateral grid surfaces (lines) divide the subsurface volume into grid 
cells. Now referring to numeral 44 of figure 5, the geologic properties of the 
subsurface volume are represented by assigning to each grid cell a particular value for 

20 each geologic property. Now referring to numeral 45 of figure 5, selected cells are 
merged to reduce the total number of cells. Now referring to numeral 46 in figure 5, 
the geologic property values are computed for the merged cells. The steps and 
various embodiments are discussed in greater detail below. 

[0053] Now referring to numeral 41 in figure 5, the input geologic data is 

25 obtained, which guides the sedimentary process simulation. The input data includes 
initial conditions, boundary conditions, and sediment discharge history. The initial 
conditions are composed of the initial (paleo) topography, depositional environments 
(deep or shallow water), and the erodability of the initial topography. The boundary 
conditions consist of inlet locations, inlet geometry, and the basin boundary flow 
30 conditions (no-flow or out-flow). The sediment discharge history consists of the flow 
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velocity, sediment-laden flow layer thickness, sediment concentration, and the 
sediment grain-size distribution function. The input data also includes information 
about locations of stratigraphic surfaces in the subsurface volume as described in 
numeral 21 of figure 2. 

« 

5 [0054] For the reference numerals 42, 43 and 44, these steps may be similar to 

the steps discussed in figure 2. In particular, with reference numeral 42, the vertical 
grid surfaces are selected in a manner that is similar to the discussion of numeral 22 of 
figure 2. For reference numeral 43, the lateral grid surfaces are selected in a manner 
that is similar to the discussion of numeral 23 of figure 2. For reference numeral 44, 
10 the geologic properties are assigned to each grid cell is similar to the discussion of 
numeral 24 of figure 2. 

[0055] Now referring to numeral 45 of figure 5, selected cells that share a 

common lateral grid surface (line) are merged to reduce the total number of cells in 
the model and thereby increase the speed of analyses that may subsequently be 
15 performed on the model. Another reason of merging the cells is to eliminate very thin 
cells that may cause numerical problems in reservoir performance simulation when 
using the geologic model. There are various criteria that can be used to determine 
which cells to merge. These criteria are clearly interrelated, and a combination of 
them maybe employed. 

20 [0056] A first approach to merging the cells may include deposit property 

changes. The deposit property changes include determining the degree of change in 
the geologic properties across the cell boundary to determine whether the cells may be 
merged. Typically, a threshold of acceptable change is selected, and cells are merged 
for cells that have changes less than the threshold value. 

25 [0057] A second approach to merging the cells may include a swept efficiency 

error (a combination effect of permeability, porosity, and other rock properties). A 
layer optimization technique can be used to determine the optimal layer-grouping 
schemes that minimize the errors of swept efficiency between the merged and original 
layers (cells). Layer optimization is a mature technology that is known to persons 
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skilled in the art and can be found in a number of publications. Examples of 
publications that disclose layer optimization include Stern, D. and Dawson, A.G.: "A 
Technique for Generating Reservoir Simulation Grids to Preserve Geologic 
Heterogeneity/ 1 paper SPE 51942 presented at the 1999 SPE Reservoir Simulation 
5 Symposium held in Houston, Texas, 14-17 February, 1999; and Li, D. and Beckner, 
B.: "Optimal Uplayering for Scaleup of Multimillion-Cell Geologic Models," paper 
SPE 62927 presented at the 2000 SPE Annual Technical Conference and Exhibition 
held in Dallas, Texas, 1-4 October, 2000. 

[0058] A third approach to merging the cells may include minimum thickness 

10 criteria. Minimum thickness criteria utilizes the thickness of the layers between 
lateral surfaces. Typically, a wide range in the thickness of the layers between lateral 
surfaces. Thin layers are generated near the edges of sand body lobes or where fine- 
grain sediments are deposited, with a sedimentary process simulation. Thick layers 
are generated near the lobe/channel centers or where coarse-grain sediments are 

15 deposited. Particularly very thin cells may cause numerical problems in reservoir 
performance simulation. A minimum thickness threshold is typically used to truncate 
very thin cells. If a cell is judged to be too thin, it may be assigned zero thickness and 
its actual thickness added to the cell above or the cell below, depending on which cell 
is closer in geological properties to the thin cell. The process is preferably continued 

20 until the thickness of the cells either exceeds the given threshold thickness or is equal 
to about zero. After merging the neighboring layers, some of layers may have zero 
thickness in some regions. Most reservoir performance simulators automatically treat 
the zero thickness cells as inactive cells. For example, Eclipse, a product of 
Schlumberger, is an example of a commercially available reservoir performance 

25 simulator that is known to persons skilled in the art. 

[0059] A fourth approach to merging the cells may be stratigraphic hierarchy. 

Stratigraphic hierarchy includes associating the lateral surfaces between cells with 
particular levels in the stratigraphic hierarchy familiar to stratigraphers. This 
hierarchy uses classifications such as "lamina", "laminaset", "bed", and "bedsef \ 
30 Cells that share surfaces below a certain level in this hierarchy may be merged For 
example, bed-scale and higher surfaces may be retained, while lamina and laminaset 
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surfaces are removed. This hierarchy and the criteria for determining which level 

♦ 

applies to a particular surface are described in Lamina, Laminaset, Bed and Bedset. 
See Campbell, C. V., Lamina, Laminaset, Bed and Bedset, Sedimentology, vol. 8 at 7- 
26 (1967). A lamina is the smallest megascopic layer that forms very rapidly from 
5 minutes to hours. A laminaset is a relative conformable succession of genetically 
related laminae bounded by surfaces of erosion, non-deposition or their correlative 
conformities, which forms rapidly from minutes to days. A bed is a relative 
conformable succession of genetically related laminae and laminasets bounded by 
surfaces (referred to as bedding surfaces) of erosion, non-deposition or their 

10 correlative conformities, which forms rapidly from minutes to years. A bedset is a 
group or set of conformable beds that differ in composition, texture, or sedimentary 
structure from those composing the beds above or below the bedset. Bedsets form 
over a longer period of time than beds. The stratigraphic hierarchy levels are 
associated with the corresponding levels of the changes of flow conditions. Avulsions 

15 and/or new inlet locations may create a high level of the hierarchy. Progradation, 
aggradation, and retrogradation may result in an intermediate level of the hierarchy 
and waves may generate a low level of the hierarchy. The stratigraphic hierarchy 
levels may be determined from the levels of the changes in the simulation conditions 
or from the inspection of the vertical distribution of the composition of the simulated 

20 deposit based on geologist's experience, as is conventionally done. The hierarchy 
levels assigned can be adjusted to achieve a favorable correlation with the cores, well 
logs, and seismic data. 

[0060] In one embodiment, once it is established that two cells may be 

merged, one of the two cells is assigned "zero" thickness and the other cell is assigned 
25 the combined thickness of the two cells. If a lateral surface contacts at least one cell 
of zero thickness at every mapview location, then it can be removed. When merging 
two cells, the choice of which cell gets zero thickness is guided by the goal of 
removing as many lateral surfaces or regions as possible. 

[0061] Figure 6(a) illustrates a hierarchy arrangement of four levels or levels 

30 1, 2, 3, and 4 labeled as 61, 62, 63, and 64 respectfully in the hierarchy arrangement 
for Hydro-Grids. Figure 6(b) illustrates the hierarchy arrangement of layers 65 within 



WO 2006/007466 PCT/US2005/021991 

-17- 

levels 61, 62, 63 and 64. The hierarchy arrangement of geologic layers is an efficient, 
need-based data structure. For example, in the exploration and early development 
stages, available data is very limited and uncertainty is high. The high-level (coarse) 
layers, such as, levels 61 or 62, of the hierarchy arrangement can be directly used to 

5 build geologic models and performance simulation models in a more efficient manner 
than the conventional geologic modeling approaches. If studies on detailed oil 
recovery efficiency are required for production decisions, the low-level (fine) layers, 
such as, levels 63 or 64, can be efficiently extracted to construct high-resolution 
models. After the desired merging process has been completed, the remaining lateral 

10 surfaces are recorded as the Hydro-Grids vertical coordinates. In some situations, the 
layer merging is not advisable because the lithology of the two layers is quite 
different. Then, the lateral grid surfaces between the two layers may be retained as a 
Hydro-Grids vertical coordinates. 

[0062] Now referring to numeral 46 of figure 5, the values of the geologic 

15 properties in every merged cell are recalculated. This determination of values for 
merged model cells is called "upscaling". Upscaling of permeability and porosity 
values is a mature technology that is known to persons skilled in the art and can be 
found in a number of publications. Examples of publications that disclose upscaling 
include Durlofsky, L. J.: "Numerical Calculation of Equivalent Grid Block 
20 Permeability Tensors for Heterogeneous Porous Media," Water Resources Research 
(May 1991) 27, No. 5, 699; Pickup, G. E., Ringerose, P. S., Forester, M. M., Jensen, 
J. L., and Sorbie, K. S.: "Permeability Tensors for Sedimentary Structures," 
Mathematical Geology (1994) 26, 227; and Li, D., Beckner, B., and Kumar, A.: "A 
New Efficient Average Technique for Scaleup of Multimillion-Cell Geologic 
25 Models," SPE Reservoir Evaluation & Engineering (August 2001), 297. 

Example 



[0063] ' Figure 7 shows a cross-section of a complex avulsive sand body 51 
composed of several sand body lobes 53 and channels 55 simulated by using a 
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physics-based sedimentary process simulator. Figure 7 shows that in complex sand 
body 51 lobes 53 and channels 55 are maintained accurately in the Hydro-Grids 58 
utilizing lateral grid lines or surfaces 54 in connection with vertical grid line or 
surface 56. The erosion (or truncation surfaces) 57 that are important to reservoir 
5 performance simulation are preserved. The color-scale 59 in Figure 7 represents 
average grain size in the phi scale (from 5 to -3). 

[0064] The Hydro-Grids idea can also be used in many commercially available 
sedimentary process simulators. Figures 8(a) and 8(b) illustrate example simulation 

10 results from a sedimentary process simulator, which populates sand bodies 61 of 
different size and shape to form a complex clastic deposit 63. Figure 8(a) illustrates 
sand body locations 61 and sizes within an example fan delta 65, and multiple lateral 
deposit time surfaces 67. Figure 8(b) illustrates multiple lateral deposit time surfaces 
67 for representative longitudinal and transverse cross-sections for sand bodies 61 in a 

15 fan delta 65. 



* 



WO 2006/007466 PCTYUS2005/021991 

-19- 



We claim: 



1. A method for geologic modeling of a sedimentary basin or reservoir, 
comprising: 

5 a obtaining geologic data; 

b. creating horizontal grid lines wherein at least one horizontal grid line 
represents deposit time lines based on the deposition rates from the geologic data; 

c. creating vertical grid lines; 

d. modeling the sedimentary properties of the sedimentary basin or 
10 reservoir from the grids of geological data created by the vertical and horizontal grid 

lines. 

2. The method of claim 1 wherein the geologic data is chosen from the group 
consisting of seismic data, outcrop study data, well log data, core data, forward and 
inverse modeling data, correlation of seismic data to sedimentary data. 

15 3. The method of claim 1 wherein the vertical grid lines are chosen to represent 
favorable vertical bins. 

4. The method of claim 1 wherein the sedimentary properties of the basin or 
reservoir are selected from the group consisting of grain size distribution, 
connectivity, net-to-gross, porosity, permeability, pore pressure and any combination 

20 thereof. 

5. The method of claim 1 wherein at the least one original horizontal grid line is 
created by estimating the deposition rate without erosion and the erosional effects are 
reconciled by adjusting at least one horizontal grid line to match the sedimentary data 



6. The method of claim 1 further comprising adjusting at least one horizontal 
25 grid line to favorably correspond to seismic reflection horizons in the seismic data. 
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7. The method of claim 6 wherein the time for adjusting at least one horizontal 
grid line to favorably correspond to a seismic reflection horizon in the seismic data 
lines is selected from the group comprising before the simulation, during the 
simulation, and after the simulation, and any combination thereof. 

5 8. A method for geologic modeling of a sedimentary basin or reservoir, 
comprising: 

■ 

a. obtaining geologic data, 

b. creating horizontal grid lines wherein at least one horizontal grid line is 
based on a deposition rate estimated from the geologic data; 

10 c. creating vertical grid lines; 

d. determining stratigraphic hierarchy levels of the horizontal grid lines, 

e. merging at least one neighboring horizontal grid line, 

f. upscaling the merged horizontal grid layers based on the stratigraphic 
hierarchy, 

1 5 g. estimating sedimentary properties from the upscaled merged layers. 

9. The method of claim 8 further comprising repeating steps (e) and (f) until all 
desired horizontal grid lines are merged. 

10. The method of claim 8 wherein the hierarchy levels assigned to the horizontal 
grid lines are selected from the group consisting of changes in the inlet conditions, 

20 changes in the simulation conditions, variation in the vertical distribution of the 
composition of the simulated deposit and any combination thereof. 



1 1 . The method of claim 8 wherein the geologic data is chosen from the group 
consisting of seismic data, outcrop study data, well log data, core data, forward and 
inverse modeling data, correlation of seismic data to sedimentary data. 
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12. The method of claim 8 herein the vertical grid lines are chosen to represent 
favorable vertical bins. 

13. The method of claim 8 wherein the sedimentary properties of the basin or 
reservoir are selected from the group consisting of grain size distribution, 

5 connectivity, net-to-gross, porosity, permeability, pore pressure and any combination 
thereof. 

14. The method of claim 8 wherein at least one original horizontal grid line is 
created by estimating the deposition rate without erosion and the erosional effects are 
reconciled by adjusting at least one horizontal grid lines to match the geologic data. 

10 15. The method of claim 8 wherein the geologic data is seismic data and further 
comprising adjusting the hierarchy levels to achieve a favorable correlation with the 
seismic data. 

16. The method of claim 15 further adjusting at least one horizontal grid line to 
favorably correspond to seismic reflection horizons in the seismic data. 

15 17. The method of claim 15 wherein the time for adjusting the horizontal grid line 
to favorably correspond to seismic reflection horizons in the seismic data lines is 
selected from the group comprising before the simulation, during the simulation, and 
after the simulation, and any combination thereof. 

18 A method for creating a geologic model of a subsurface volume, comprising: 
20 a. obtaining geologic data; 



b. 



creating vertical grid surfaces; 



c. 



creating lateral grid surfaces that correspond to surfaces of geological 
time during the deposition of sediments in a subsurface volume; and 



25 



d. 



representing geologic properties of the subsurface volume as values in 
each cell created by intersections of the vertical gird surfaces and the 
lateral grid surfaces. 
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19. The method of claim 18 wherein the geologic data comprises at least one of 
seismic data, outcrop study data, well log data, core data, numerical modeling data, 
stratigraphic surfaces interpreted from seismic data, and any combination thereof. 

20. The method of claim 18 wherein the vertical grid surfaces partition the 
5 subsurface volume into regions based on lateral variation of geologic properties. 

21. The method of claim 20 wherein the partitioning of the regions is configured 
to have large bins for little lateral variation in geologic properties and small bins for 
large lateral variation in geologic properties are divided into small bins. 

22. The method of claim 18 wherein the geologic properties of the subsurface 
10 volume comprises at least one of grain size distribution, connectivity, net-to-gross, 

porosity, permeability, pore pressure and any combination thereof. 

23. The method of claim 18 wherein at least one of the lateral grid surfaces is 
adjusted to match the geologic data. 

24. The method of claim 18 wherein at least one of the lateral grid surfaces is 
15 adjusted to account for the removal of sediment by erosion subsequent to deposition. 

25. The method of claim 18 wherein the lateral grid surfaces are determined by 
numerical simulation of sedimentary processes based on the geologic data. 

26. A method for creating a geologic model of a subsurface volume, comprising: 
a. obtaining geologic data; 

20 b. creating vertical grid surfaces; 

c. creating lateral grid surfaces corresponding to surfaces of geological 
time during deposition of sediments in a subsurface volume; 



25 



d. representing geologic properties of the subsurface volume as values in 
each cell created by the intersections of the vertical grid surfaces and the 
lateral grid surfaces; 
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e. merging selected grid cells; and 

f. calculating geologic property values in the merged grid cells. 

27. The method of claim 26 further comprising repeating steps (e) and (f) until 
each of the desired grid cells are merged. 

5 28. The method of claim 26 wherein grid cells are selected for merging based on 

< 

criteria including at least one of a classification of a common grid surface in a 
stratigraphic hierarchy system, grid cells having substantially equal geologic property 
values, one of the grid cells having a thickness less than a threshold value, the grid 
cells merging with minimum swept efficiency errors, a flow properties during their 
10 deposition having been about equivalent 

29. The method of claim 26 wherein the geologic data comprises at least one of 
seismic data, outcrop study data, well log data, core data, numerical modeling data, 
stratigraphic surfaces interpreted from seismic data, and any combination thereof. 

30. The method of claim 26 wherein the vertical grid surfaces partition the 
1 5 subsurface volume into regions based on geologic properties of layers. 

31. The method of claim 30 wherein the partitioning of the regions is configured 
to have large bins for little lateral variation in geologic properties and small bins for 
large lateral variation in geologic properties are divided into small bins. 

32. The method of claim 26 wherein sedimentary properties of the subsurface 
20 volume are selected from at least one of grain size distribution, connectivity, net-to- 
gross, porosity, permeability, pore pressure and any combination thereof. 

33. The method of claim 26 wherein at least one of the lateral grid surfaces is 
adjusted to match the geologic data. 

34. The method of claim 33 further adjusting at least one of the lateral grid 
25 surfaces to favorably correspond to seismic reflection horizons. 
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35. The method of claim 26 wherein the lateral grid surfaces are determined by 
numerical simulation of sedimentary processes based on the geologic data. 

36. The method of claim 26 wherein at least one of the lateral grid surfaces is 
adjusted to account for removal of sediment by erosion subsequent to deposition. 
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